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Abstract 


In  this  paper,  we  consider;,  the  linear  quadratic  optimal  control 
problem  on  infinite  time  interval  for  linear  time-invariant  systems  defined  on 
Hilbert  spaces.  The  optimal  control  is  given  by  a  feedback  form  in  terms  of 
solution  n  to  the  associated  algebraic  Riccati  equation  (ARE).  A  Ritz  type 
approximation  is  used  to  obtain  a  sequence  r$  of  finite  dimensional 
approximations  of  the  solution  to  ARE.  A  sufficient  condition  that  shows  nN 

converges  strongly  to  n  is  obtained.  l/nder  this  condition,  we  derive  a 

1 

formula  which  can  be  used  to  obtain  a  rate  of  convergence  of  nN  to  n. 
We  demonstrate  and  apply  the  results  for  the  Galerkin  approximation  for 
parabolic  systems  and  the  averaging  approximation  for  hereditary  differential 
systems. 


Assume  Z,  U  and  Y  are  Hilbert  spaces.  Consider  the  evolution 


equation  on  Z 

(1.1)  z(t)  =  A  z(t)  +  B  u(t)  ,  z(0)  -  z0  €  Z 

where  u(t)  is  T  alued  control  function,  A  is  the  infinitesimal  generator  of 
strongly  continue  emigroup  S(t)  on  Z,  and  B  €  2(U,Z).  The  Y-valued 
observation  function  y  is  given  by 

(1.2)  y(t)  =  C  z(t)  ,  t  >  0  . 

We  assume  that  C  €  t( Z,Y).  We  interpret  the  equation  (1.1)  in  the  mild  sense; 
the  solution  of  (1.1)  is  given  by 

(1.3)  z(t)  «  S(t)zQ  +  f  S(t-s)R  u(s)ds  . 

Consider  the  minimization  problem;  minimize  the  cost  functional 

(1.4)  J(u,z0)  «  [  (||y(t)||J  +  j|u(t)||J)dt 

•’o 

subject  to  (1.3).  Then  the  following  result  is  well-known  [10],[11]: 

Theorem  1.1  Assume  ( A.B )  is  stabilizable  and  (A.C)  is  detectable.  Then 

there  exists  a  unique  nonnegative  self-adjoint  solution  n  to  the  algebraic 
Riccati  equation  in  Z: 

(1.5)  (/l*n  +  ru  -  nss*n  +  C*C)z  *  0  for  all  z  €  dom(A ) 
and  the  optimal  solution  u°  to  (1.4)  is  given  by 

u°(t)  =  -5*n  T(t)z0 

where  T(t)  is  the  strongly  continuous  semigroup  generated  by  A  -  BB*n 
and  it  is  uniformly  exponentially  stable. 


Here  we  have 


Definition  1.2  (1)  ( A,B )  is  stabilizable  if  there  exists  an  operator 

K  €  E(Z,U)  such  that  A-BK  generates  a  uniformly  exponentially  stable 
semigroup  on  Z. 

(2)  ( A.C)  is  detectable  if  there  exists  an  operator  G  €  £(Y,Z)  such  that  A  -  GC 
generates  a  uniformly  exponentially  stable  semigroup. 


The  purpose  here  is  to  construct  a  finite  dimensional  approximation 
of  the  optimal  feedback  gain  operator  £?*n.  Let  us  consider  a  sequence  of 
approximating  problems  (ZN,/4N,fiN,CN),  let  ZN  be  a  sequence  of  finite 
dimensional  subspaces  of  Z  and  PN  be  the  orthogonal  projection  of  Z  onto 
ZN.  Assume  /1N  :  ZN  -ZN  ,  Z?N  :  U  -*  ZN  and  CN  :  ZN  -  Y  are  continuous. 
Then  consider  the  Nth  approximating  problem  of  (1.4) 


(1.6)  minimize  J 


subject  to 


N(u,z0)  =  J“(|jCNzN(t)||2  +  |ju(t)||2)dt 


(1.7)  zN(t)  =  SN(t)PNz0  +  [  SN(t-s)BNu( 

*  n 


/here  SN(t)  =  cA  t  >  0.  Then  the  optimal  control  uN  of  (1.6)  is  given 


uN(t)  =  .gN*nNeHN-5NfiNnN)tpNZo  f  t  >  0 

where  nN  :  ZN  -•  ZN  is  self-adjoint  and  satisfies  the  Nth  approximating 
algebraic  Riccati  equation  in  ZN; 

(1.8)  /1NV  +  n n/4n  -  nNRN5N*nN  +  cN*cN  «  o  . 

Here,  fiN*nN,  n  >  1  yields  a  sequence  of  finite  dimensional  approximations  of 


//VvvVV/v'.v-/.'.-  v.'-w-v  av.v.v  .v.v;.-, 


the  optimal  feedback  gain  [3]. 

In  this  paper  we  first  obtain  a  condition  on  (ZN,/4N,flN,CN)  for  which 
(1.8)  admits  a  unique  nonnegative  solution  nN,  and  nNPN  converges  strongly  to  n 
in  §2.  Such  a  condition  has  been  discussed  in  [2],  [3]  but  the  condition  in  this 
paper  improves  those  in  [2],  [3],  i.e.,  we  introduce  the  uniform  detectability 
condition  (see,  (H3)  in  §2,  for  the  definition)  which  is  additional  to  those 
considered  in  [2],  and  using  this  condition,  we  are  able  to  show  that  there 
exists  an  integer  N0  such  that  for  N  *  N0 

|| e(^N-BN5N  nN)tpN||  $  m  e*u  1  ,  t  >0 

for  positive  constants  M  >  1  and  u  (independent  of  N  )  N0).  This  assertion 
is  a  part  of  assumptions  in  [2,  Theorem  2.2],  The  uniform  detectability 
condition  is  satisfied  if  C*C  is  coercive,  which  is  assumed  in  the  discussions  in 
[2,p.  693],  Thus,  the  uniform  detectability  condition  can  be  regarded  as  a 
relaxation  of  the  coercivity  assumption  mentioned  above.  Next,  under  the 
condition  in  §2  we  derive  a  formula  which  provides  a  rate  of  convergence  of 
nN  to  n  and  apply  the  formula  for  specific  examples. 


We  assume  the  following.  Lei  5N(t)  =  e^1 
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CUl 


,  t  >  0 


(HI)  For  each  z  €  Z,  we  have 

(i)  SN(t)PNz  -  S(t)z  ,  and 

(ii)  SN(t)*PNz  -  S*(t)z  , 

where  the  convergences  are  uniform  in  t  on  bounded  subsets  of  [0,”°). 

(H2)  (i)  For  each  u  e  U,  fiNu  -*  B u  and  for  each  z  c  Z 

flN*PNz  -  B*z 

(ii)  For  each  z  €  Z,  CNPNz  -*  Cz  and  for  each  y  t  Y 
CN*y  -  C*y. 

(H3)  (i)  The  family  of  the  pairs  (/4N,BN)  is  uniformly  stabilizable:  i.e. 

there  exists  a  sequence  of  operators  KN  c  I(ZN,U)  such  that 


sup||KN||  <  ®  and 


||e(^N-fiNA:N)tpNjj  <  Mje'^i1  ,  t  >  0 

for  some  positive  constants  Mj  >  1  and  Uj  . 

(ii)  The  family  of  the  pairs  (/4N,CN)  is  uniformly  detectable;  i.e. 
there  exists  a  sequence  of  operators  CN  c  £(Y,ZN)  such  that 
sup||CN||  <  ®  and 

|jC(^N-CNCN)tpN||  <  M2e"^jt  ,  t  >  0 
for  some  positive  constants  M2  >  1  and  u>2. 

(1)  Suppose  fiN  =  PNfl  and  CN  -  CPN.  Then  (H2)  holds 


Remark  (1) 


since  it  follows  from  (HI)  that  PNz  -•  z  for  all  z  t  Z. 

(2)  The  assumption  (H3)  is  closely  related  to  the  preservation  of 

exponential  stability  under  approximation  in  [3, Conjecture  7.1]  and  it  is  shown 
in  [2]  that  (H3)  (i)  ((POES)  in  (2])  is  satisfied  for  parabolic  systems  using  the 
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Galerkin  approximation. 

(3)  A  natural  way  to  argue  (H3)  is  to  take  K N  =  K PN  and  CN  =  PNC  for 
some  K  £  £(Z ,U)  and  G  t  £(Y,Z)  such  that  A  -  BK  and  A  -  GC  generate 
uniformly  exponentially  stable  semigroups  on  Z. 

Theorem  2.1  Suppose  (H1)-(H3)  are  satisfied.  Then  for  each  N,  (1.8) 

admits  a  unique  nonnegative  solution  nN,  sup||nN||  <  and  there  exist 
positive  constants  M,  >  1  and  u,  (independent  of  n)  such  that 


e(^N-fiNfiN  nN)tPNl|  <  M,e'W»l  ,  t  >  0 


Proof:  The  proof  is  based  on  the  arguments  in  [II],  The  existence  and 
uniqueness  of  solutions  to  (1.8)  follow  from  Theorem  1.1.  Since 
<nNPNz,z>  =  min  JN(u,z),  (H3)  (i)  implies  that 


<nNPNz,z  >  <  JN(-tfNzN(  ■  );z) 


f"  (||CNe^N*BN^N>tPNz||J  +  ||^Ne^N-5N^N)tPNz||2)dt 

*  n 


<  Bllzlt2  for  some  positive  constant  B 


Since  nN  is  self-adjoint,  nonnegative  definite,  this  implies  that  ||nN||  4  & 
By  the  variation  of  constants  formula 

(2.1)  C(/1N  -fiNSN  nN)t  =  TN(t)  +  f  xN(t-s)(C7NCN-5N5N*nN)e(/,N'flNfiN  nN)s  ds 

•  n 


where  TN(t)  =  c(^N-CNCN)t>  t  >  0  Here,  from  (1  8) 

MN-f?NBNV)*nN  +  nNMN-sNBNV)  +  nN£Nz?NV  +  cN*cN  =  o  , 

so  that  if  zN(t)  «  e(AN-BNBN  rfJ)tpNZ)  t  >  q,  then 
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1  <zN(t),nNzN(t)>  +  ||SN*nNzN(t)||J  +  j|CNzN(t)  ||J  =  o 
Thus,  for  all  t  >  0 

(2.2)  <nNzN(t),zN(t)>  +  [  (  )]SN*nNzN(t)j|2  +  ||CNzN(t)||2)dt 

<  <nNPNz,z>  <  6||z|[2 

Now,  from  (2.1),  we  have  for  all  t  >  0 

f  IN’  +  P  (|pN||!  *  |*N|f>  |t(||sN*nNzN(s)||I 

J0  W2  Uj  J  0 

+  ||CNzN(s)||2)ds 

where  we  have  used  the  Young’s  inequality.  From  (2.2),  we  have 

r  m2 

||zN(t)||2dt  <  — L(u2  +  26  Sup(||CN||2  +  ||5Nj|2))||z||2 
*  0  u2 

for  all  z  €  Z.  Therefore,  the  theorem  follows  from  the  Datko’s  theorem  [7], 

(Q.E.D.) 

The  following  is  a  consequence  of  [3,  Theorem  6.9]  and  [2,  Theorem 

2.2], 

Corollary  2.2  Suppose  (A,B)  is  stabilizable  and  ( A,C )  is  detectable  and 

assume  (HI)  -  (H3)  are  satisfied.  Then  the  unique  nonnegative  solution  nN 

to  (1.8)  converges  strongly  to  n,  the  unique  solution  to  (1.5). 

Theorem  2.3  Suppose  that  B  is  compact  and  Bs  =  P NB  and  that 

(HI)(i)  and  (H3)(i)  are  satisfied.  Then  (A,B)  is  stabilizable. 

Proof :  Let  us  consider  the  case  C  *  I  and  CN  =  PN  with  Y  =  Z.  Then  it  is 
easy  to  show  that  ( A,C )  is  detectable  and  (^N,CN),  N  >  1  are  uniformly 
detectable  since  (Hl)(i)  implies  that  for  some  M  ?  1  and  u  independent  of  N, 


t &*** 


||SN(t)PN||  «  Me^*,  t  )  0.  It  then  follows  from  Theorem  1.1  and  (H3)(i)  that  for 
each  N,  (1.8)  with  CN  =  PN  has  a  unique  solution  fiN.  Using  the  same 
argument  as  in  the  proof  of  Theorem  2.1,  we  have  ||n  )|  i  B  for  some  positive 
constant  B.  Thus,  by  Theorem  6.5  in  [3],  there  exists  a  subsequence  of  nN 

A 

converges  weakly  to  some  nonnegative,  self-adjoint  operator  n  We  will  show 
that  n  satisfies  (1.5)  with  C  =  I.  Since  PNnN  =  nN,  BN*nN 

s*pNnN  =  b*  n*  Since  B*  is  compact,  for  each  z  e  Z,  B  ■>  ri  jPNz  converges 

A  A 

strongly  to  S*nz.  It  now  follows  from  [3,  Theorem  6.7]  that  n  satisfies  (1.5) 

A 

But  since  ( A.C )  is  detectable,  by  [10,  Theorem  3.2],  A  -  BB* n  generates  a 
uniformly  exponentially  stable  semigroup  on  Z. 


(Q.E.D.) 

Remark  2.4  Roughly  speaking,  Theorem  2.3  means  that  the  uniform 

stabilizabil  ity  implies  the  stabilizability  of  ( A.B ).  The  dual  statement  of 
Theorem  2.3  also  holds:  i.e.,  suppose  C  is  compact,  CN  =  CPN,  then  (Hl)(iij 
and  (H3)(ii)  imply  that  (A.C)  is  detectable.  This  statement  can  be  proved  by 
applying  the  exactly  same  arguments  as  in  the  proof  of  Theorem  2.3  to  the 
dual  Riccati  equation 

(.41  +  LA*  -  LC*Cl  +  I)  z  =  0 
for  all  z  €  dom(/l*). 


i  -A 


In  this  section,  we  assume  that  (HI)  and  (H3)  are  satisfied  and  let 


SN  =  pN5  and  CN  =  CPN.  Moreover,  we  assume 
(H4)  For  each  z  «  Z,  nz  e  dom(A*)  and  B  is  compact. 
From  (1.5),  we  have  for  all  z  e  dom (A) 


2<nz,-4z>  -  <fl*nz,fl*nz>  +  <Cz,Cz>  =  0 
Thus,  (H4)  and  the  density  of  dom(/l)  in  Z  imply  that  for  all  z  t  Z 


2</4*nz,z>  -  <B*Ilz,B*nz>  +  <Cz,Cz>  =  0  . 

Define  the  self-adjoint  operator  nN  =  PNnPN.  Then  for  all  x  c  ZN 


(3.1)  2<04N*nNx,x>  -  <B*nNx,  B*nNx>  +  <CNx,CNx>  +  <Anx,x>  =  0  , 


where  an  e  I(ZN)  is  a  self-adjoint  operator  defined  by 


(3.2)  <Anx,x>  =  2 <(A*  -  /iN*PN)nx,x> 


+  <fi*(nN  -  nN)x,S*(riN  +  nN)x>  for  all  X  e  ZN 


From  (1.8),  for  x  €  Zr 


(3.3)  2<C4N*nNx,x>  -  <5N*nNx,£N*nNx  >  +  <CNx,CNx>  =  0 


Hence  by  subtracting  (3.1)  from  (3.3) 


2  <mn  -  fiNsN*nN)x,  (nN  -  nN)x> 


+  <sN*(nN  -  nN)x,  flN*(nN  -  nN)x>  -  <anx,x>  -  o 


for  all  x  «  ZN.  Or  equivalently 


(3.4)  nN  -  nN  =  f  eMN-sNsNV)t((nN  .  n")B*B"\n*  -  it*)  -  an) 


°MN-BNsNV)t  dt 


Similarly,  subtracting  (3.3)  from  (3.1),  we  obtain 


(3.5)  nN  -  nN  =  f  cMN-sNaN*nN)t((nN  .  .  ft*)  +  an} 

J  n 


X  B  nN)tdt  . 

Here,  from  Theorem  2.1,  we  have 


nN>t  pN  ji  <  M  .-v,  t  >  o 


u.t 


with  Mj  >  1  and  ws  >  0.  Since  nNPN  -*  n,  strongly  by  Corollary  2.2  and 
B  is  compact 


sN  (nN  -  nN) |[  *  jj(nN  -  n)PN5j|  -  0  as  n  - 


N  .  nlDNj 


Hence  by  the  variation  of  constants  formula  and  the  Gronwall’s  lemma, 

pN  (!  <  Mje(*us+  llfi  II  ||aN*(nN-nN)  ||)t 

It  then  follows  that  there  exists  an  integer  N0  such  that  if  N  >  N0, 


.(^N-flNfiN*riN)t  pN 


||  <  M,e  2  ,  t  >  0  . 


Now,  from  (3.4)  for  all  x  t  ZN 


(3.6) 

<(nN  -  nN)x,x>  -  J 

j|flN(nN  -  fiN)eMN-aNaNV)tX|j2  dt 

0 

-  ■ 

.  f"  <c(/<N-sNfiN*nN)tx^  ANcuN-BNfiN*nN)tx>dt 

and 

from  (3.5) 

*  0 

(3.7) 

<(nN  -  nN)x,x>  -  j 

jjSN(nN  *  nN)e(/*N'flNfiN  "N>lx |(*dt 

0 

- 

f  <cUN-BNBN*nN)tx  ANC(^N.BNBN-N)tx>dt 
^0 

These  inequalities  imply  that  for  x  «  ZN 
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,  <(nN  -  nN)x,x>|  <  — 


N"»x||2 


so  that 


k/--‘ 


y.  ■' 


(3.8)  ||nN  •  fiN||  <  |A"| 

s 


where  from  (3.2) 


(3.9)  ||An||  <  2 1| (A*  -  >4N*PN)nj|  +  20||fl||||(B*-SN*)nj|  for  all  N  »  N„. 


In  this  section  we  discuss  the  examples  in  which  (H1)-(H4)  are 
satisfied  and  then  apply  the  formula  (3.8)  and  (3.9)  to  obtain  a  convergence 
rate  of  nN  to  n. 


Assume  V  and  H  are  Hilbert  spaces  and  V  C  H  with 

continuous  dense  injection  i.  Consider  a  bilinear  form  o  on  V  such  that 


(4.1)  |o(u,v)|  i  cllull  jlvll  for  u,v  e  V 


(4.2)  o(u,u)  >  u||u||2-  pj|u||2  for  u  e  V 


where  u  >  0.  It  then  follows  from  [9]  that  there  exists  an  operator  A  e 
r(V,V*)  such  that 


(4.3)  o(u,v)  =  <-^u,v>v#v  for  u,v  e  V  , 

where  V  c  H  =  H*  C  V*  and  H  being  the  pivoting  space,  and  that  A 

on  H  with 


(4.4)  dom(/4)  =  (x  €  H  :  Ax  e  H}  dense  in  V, 

generates  the  analytic  semigroup  on  H  and  V*.  For  given  B  e  S(U,H)  and 
C  «  2(H,V)  consider  approximating  problems  (ZN,/4N,BN,CN);  i.e.  let  ZN  be  a 
sequence  of  finite  dimensional  subspace  of  V  and  AN:  ZN  -*  ZN  is  defined  by 

(4.5)  <-/4nz,x>  =  o( z,x)  for  z,x  €  ZN  . 

Let  PN  be  the  orthogonal  projection  of  H  onto  ZN  and  assume  BN  = 
PNS  and  CN  =  CPN.  We  assume  the  approximation  condition: 


-12- 


For  each  z  £  V  there  exists  an  element  zN  «  ZN 
(4.6)  such  that  llz  -  zN II  <  «(N)  where  €(N)  -*  0  as  N  -  00  . 


It  then  follows  from  [2]  that  (HI)  holds  and  if  ( A.B )  is  stabilizable  and  (A,C)  is 
detectable,  then  (H3)  holds.  Thus  from  Corollary  2.2,  nN  converges  strongly  to 
IL  However  one  cannot  apply  the  formula  (3.8)-(3.9)  as  it  is,  since 
nZ  C  dom(/l*)  is  the  maximal  regularity  without  assuming  any  regularity  of  C. 
This  can  be  demonstrated  by  the  following  example.  Consider  the  case  when 
H  =  L2(0,1)  and  V  =  Hj(0,l),  and 


f  d  d  , 

(u,v)  =  - —  u(x)  — —  v(x)dx  for  u,v  €  H* 

Jn  dx  dx 


Let  us  consider  the  Liapunov  equation  on  H 


(4.7)  Al  +  ZA  +  Q  =  0 

where  Q  is  self-adjoint  operator  on  H.  If  for  each  z  £  Z, 


(Ez)(x) 


=  f  «x, 

•  n 


y)z(y)dy  and  (Qz)(x) 


f  Q(x, 

*  n 


y)z(y)dy. 


then  <p  satisfies  A<£  +  q  =  0  with  Dirichlet  boundary  condition, 

3 2  32  , 

where  A <£  =  - 0  #  fo'  4>  £  H  ([0,1]  x  [0,1]).  In  general  (e.g.,  see 

3x2  3y2 

[6], [8]) 

[  f  \\l_i  2  +  |_L_<J>|J]  dxdy  <  M  f  f  I q  1 2  dxdv  . 

Jo  Jo  ^3x2  'ay2  1  J  Jo  Jo'  1 

This  implies  EL2  C  dom(/4)  is  the  maximal  regularity. 

Hence  we  will  modify  the  arguments  in  Section  3  to  improve  the 

formula  (3.8)-(3.9)  for  this  example.  First  we  note  that  in  (3.2)  for  x  £  ZN 

|<M*-.4NV‘)nx,x>|  -  |o(x,(nN  -  n)x)| 

<  c||x|y|<n  -  fi")xj|v  by  (4.1) 


Thus  from  (3.2) 


I  <Anx,x  >  I  ||(n-nN)x  ||v 

(4.8)  7  =  sup  -! - L  <  2  c  sup  JJ - — 


lx  II  *0  ||x  | 


x*0  || x  | 


Mg*l  sup  Egjfelr , 


x*0 


where  ||x||h  <  «||x  ||v- 


Lemma  4.1  There  exists  a  positive  constant  M  such  that 

f*  oNoN*nN\f  m2  2 

i  He'-4  -do  n  )t  pNx  |j  <  M  |jx  |J  ,  and 
J  n  V  H 


i; 


(/N  oNdN^Nw  v  2  2 

\\',A  B  ""  pN>llv  d'  <  MIIMI„ 


Proof:  Let  tN(t)  -  e</,N-fiNsNnN>t  PNx  ,  t  >  0.  Then  lN(t)  satisfies 

~  *N(t)  *  (^N-5NJ?N*nN){N(t)  ,  t  >  0  , 
dt 

so  that  from  (4.5) 

L  JL  pN(t)||2  +  oUN,SN)  «  -  <flN*lN(t),BNVtN(t)> 

2  dt  11  "h  h 

and  from  (4.2) 

7  *  “ll'NwiC  ‘ <p  *  28»sII!)  lllN(,)»l 

The  integration  of  this  inequality  with  respect  to  t  yields 

7  ll'N(,)C  *  “  f0  ll'N(,»llv  ‘  7  l'"(0)l*  *  (p  *  2e»Bl|!)  f0  II' 


Now  the  lemma  follows  from  Theorem  2.1. 
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1 

U 

m 


I 

D 


1 


It  then  follows  from  Lemma  4.1  and  (4.8)  that 


f"  <an  e(/4N-fiNBN*nN)tx  ^  e(/lN-flNfiN*nN)tx  > 

•  n 


t~  ,jN  pN  dNVNw  w  2  2 

<  7  I|e(/4  'B  B  n  PNX  ||  dt  <  M7||x||  . 

J  n  V  ^ 


,  iN  dN pN  nNw 

Similary  for  e(/4  b  11  '  ,  t  >  0.  Therefore  we  obtain,  using  (3.6)  and 


(4.9)  ||nN  -  nN  ||  <  M7  • 


where  y  is  given  by  (4.8). 


Consider  the  (1-dimensional)  parabolic  control  system  [2]; 


0  0  0  0 

—  z(t,x)  =  _  (p(x)  —  z)  +  q(x)  —  z  +  r(x)z 
0t  ox  ox  ox 


+  I  b|(x)uj(t)  in  (0.1) 
i=i 


with  boundary  condition  z(t,0)  *  z(t,l)  -  0,  where  p  €  C^O,!),  being  bounded 


below  bv  a  positive  constant  w,  —  q,  r  €  L  (0,1),  and  b:  €  L2(0,1),  i  =  l,....m. 

dx 


In  this  case,  H  =  L2(0,1)  and  V  «  Hl(0,l),  and  the  bilinear  form  o  is 


given  by 


f1  d  d  d 

o( U,v)  =  [p(x)  — u  —  V  -  (q(x)  — u  +  r(x)u)v]dx 

Jn  dx  dx  dx 


B:  ■*  LJ(0,1)  is  defined  by 


(Bu)(x)  *  I  bj(x)uj  for  u  (  P"1  , 

i=l 


and  dom(/4)  *  do m(A*)  «  HJ(0,1)  n  Hi(0,I).  Let  us  consider  the  following 


finite  dimensional  subspace  ZN  of  V: 


ZN  =  {z  €  H  :  z(x)  -  Z  OjB^  (x)  ,  bj  c  R) 


where  BN(  • ).  i  *=  1 . N-l  are  the  linear  B-spline  elements  on  the  interval  [0,!]; 

i.e., 


i+1 

-N(x- - )  , 

N 


N(x-  ■ 


i-1 

N~ 


elsewhere 


Then  the  approximation  condition  (4.6)  is  satisfied  [8].  Suppose  (A.B)  is 

stabilizable  and  ( A.C )  is  detectable.  Then  (1.5)  has  the  unique  solution  n 

and  using  a  similar  arguments  to  those  given  above  to  show  the  regularity  of 

solutions  to  Liapanov  equation  (4.7),  one  can  show  that  for  x  «  H,  fix  € 

dom(.4*).  Since  A*  is  closed  in  H  and  dom(/!*)  C  V,  by  the  closed  graph 

theorem,  there  exists  a  positive  constant  k.,  such  that  linzll  .  t  k.||z||  ■> 

1  11  111  "L‘(0,1) 

Hence  the  fundamental  error  estimate  (e.g.,  (8J)  gives 


||TU  -  (lN2||L,  <  k, 

II"2  '  ""'Ll  ‘  k> 
H0 

for  some  positive  constants 


[n]  BzBlj 

kj.kj  Now  it  follows  from  (4.8)  and  (4.9)  that 
for  some  constant  k. 


4.2  Hereditary  Differential  Systems 

Consider  the  hereditary  differential  system  in 

x(t)  =  A0x(t)  +  AjX(t-r)  +  j  A(0)x(t+8)de  +  Bu(t) 

*  -r 

x(0)  «  n  and  x(0)  -  ^(6),  -r  <  0  <  0 


(4.10) 


and  the  optimal  control  problem;  for  given  initial  data  z  =  (n,$)  «  R"  x 
L2(-r,0;  R")  ,  minimize  the  cost  functional 

-QD 

(4.11)  J(u,z)  =  |  (|y(t)|J+  | u(t) 1 2)dt  . 

Here,  x  (  PP,  u  t  R"1  and  y  t  Ik’  and  the  element  of  A(  )  is  square 
integrable.  It  is  shown  [1]  that  (4.10)  and  (4.11)  are  equivalently  formulated  as 
the  problem  (1.1)  -  (1.4)  on  the  product  space  Z  «  fi"  *  L^-r^R");  i.e.,  z(t)  = 
(x(t),x(t  +■))  e  Z  is  the  mild  solution  of  (1.1)  with 

dom(/4)  =  {(h,$)  e  Z  :  $  €  H^-r.O)  and  $>(0)  =  b)  , 

for  (<$K0),<f>)  €  dom(4) 

<4(«0),«  =  (Ao&0)  +  A^-r)  +  j  A(6)<«e)de,  0)  . 

J  -r 

The  input  operator  B  :  R™  -•  Z  and  the  output  operator  C  :  Z  -  R*3  are 
given  by 

B U  *  (Bu,0)  €  Z  and  C(n,4>)  «  Cb  . 

Let  us  consider  the  averaging  approximation  [1]  of  (4.10);  let 

N 

ZN  =  {z  €  Z  :z=  (a0, 1  akx  s  ),  ak  €  Rn  ,  0  <  k  i  N)  c  Z  , 

k=l  [- — r, - r) 

N  N 


and  -4n  has  the  matrix  representation  (Qn)''Hn  on  p"(N+1)  when  ZN  is 
identified  with  R"^N+1^  by  its  coordinate  vector  col(aJ,...,aJ),  where  the  block 
diagonal  matrix  QN  and  the  block  Hessenberg  matrix  HN  are  given  by 


I 

r 

[an 

AN  aN  1 

qn  = 

—  I 

N 

and  Hn  - 

I 

-I 

r 

—  I  . 

I  -I 

with  Aq  =  A0  ,  Aj  =  N  A  (0)d0  and  A^J  =  At  +  n  A(0)d0.  Note  that 

J --r  -r 

N 

PNS  =  B  and  CPN  =  C.  Set  =  B  and  CN  =  C.  Then  (/1N)*  has  the 

matrix  representation  (QN)_1HN  on  Rn(N+1),  (Hl)(i)  is  proved  in  [1]  and 
(HI )(ii)  is  proved  in  [3].  Using  the  arguments  in  [5],  [7]  one  can  show  that 

(H3)  is  satisfied  (i.e.,  (i)  is  straightforward  but  (ii)  is  not  so).  Thus,  the 

formu.a  (3.8)-(3.9)  yields 

|jnN  -  nN ||  (  2|p*  -  /4N*PN)nj| 

By  the  regularity  result  in  [4],  if  A(  )  «  H1(-r,0;Rnxn),  then 
A* n  +  C*C  e  dom(/l*) 

where  dom(/l*)  =  {(y,0)  £  Z  :  0  £  H1  and  0(-r)  =  A^y}  and  A*( y,0)  = 

WO)  +  Ajy,  -  0(0)  +  At(  ) y)  c  Z  [3],  Since  C*C(n,0)  =  (CTCn,0)  €  Z  for 
(n,0)  £  Z,  this  implies  that  if  nz  =  (y,0),  then  0  €  H1  so  that  0  t  H2,  and 

since  A*  is  closed,  ||0||  <  M||z||  f°r  some  constant  M  It  then  follows 

V  z 

from  the  arguments  and  error  estimate  in  [1],  [3]  that 

|p*-/lN*PN)n(y,0)||  <  —  ( |  y  |  +  ||0||  ,) 

''N  h 

||  nN  -  nN||  =  0(  L)  . 


Hence  we  obtain 
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